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BACKGROUND 


• CMOTT group at LeRC has been in existence for about 4 years. In the 
first 3 years, its main activities were in developing and validating tur- 
bulence and combustion models for propulsion systems, in am effort to 
remove the deficiencies of the existing models. Two workshops on com- 
putational turbulence modeling were held at LeRC (1991, 1993). 

• A peer review of turbulence modeling activities at LeRC was held in 
September, 1993. Seven peers (GE, P&W, RocketDyne, Cornell, Berke- 
ley and NASA Ames) conducted the peer review. The objective of the 
peer review was to assess the turbulence program at LeRC/CMOTT 
and to suggest the future direction of turbulence modeling activities for 
propulsion systems. 

• Important messages from the peer review: 

0 “LeRC should spend substantial effort being responsive to indus- 
try’s current pressing perceived needs; this involves extensive dis- 
cussion with industry during every phase of model development, 
analysis of industry’s problems, goal oriented model development, 
evaluation of models relative to industry’s intended application 


C 1 “LeRC has an obligation not only to respond to industry’s re- 
quests for help, but to play an autonomous, independent leader- 
ship role in providing models of the highest quality, which can 
be employed not only by the aircraft gas turbine and rocket indus- 
tries but also by other industries . . ’.’ 

<> “In the present financial climate, industry does not have the re- 
sources to undertake model development and evaluation. LeRC’s 
help in this regard via the creation of its turbulence modeling ef- 
fort, is, therefore, welcome from the industry’s standpoint.” 

§ “It is important to work with the industry to evaluate the models 
and rank-order them by performance and cost in order to identify 
the most appropriate models for particular situations.” 

0 Many other useful suggestions and comments including collabora- 
tion with industry, joint programs, industry-wide workshop . . . 
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PROGRAM GOALS AT CMOTT 


• Develop reliable turbulence (including bypass transition) and combustion 
models for complex flows in propulsion systems 

• Integrate developed models into deliverable CFD tools for propulsion 
syst ems in collaboration with industry. 

PROGRAM APPROACH 


• Develop turbulence and combustion modules for industry customers 

• Industry collaboration and technology transfer 

• Model development for propulsion systems 

0 One-point moment closures for non-reacting flows 
0 Scalar PDF method for turbulent reacting flows 
0 Valida tion of existing and newly developed models 
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Development of 

Turbulence and Combustion Modules 


Objective 

0 Build a quick and efficient vehicle for technology transfer to indus- 
try 

The features of the turbulence module: 

0 It contains various turbulence models from which users can choose 
the appropriate model for flows of interest 

0 It is self-contained, i.e., it contains its own solver for turbulence 
model equations 

0 It can be easily linked to industry’s CFD codes 

Turbulence module for NPARC code has been developed, tested, and is 
ready to be released 

0 The models built-in at the present time: 

Mixing length, Chien k - e, CMOTT k - e models 
0 The model to be built-in: 

CMOTT algebraic Reynolds stress, Reynolds stress transport 
equation models and other models based on the request from in- 
dustries. 

Built-in robust, realizable numerical solver for model equations. 
General turbulence modules 

0 Can be used for both compressible and incompressible flows. 

0 Interface programs for different industry CFD codes 
Built-in models will be periodically updated. 
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H H O S O 


Collaboration with Industry 
and Technology Transfer 

• Joint research programs with industry 

0 Preliminary programs with engine companies and others have been 
initiated (GE, P&W, RocketDyne, Naval Research Laboratories) 

^ Develop further joint research programs related to the industry’s 
projects 

• Industry-wide workshops will be a regular program (once every two 
years) 

<£> Release Lewis turbulence and combustion modules to industries 
0 Discuss the needs of industry 

Models developed at CMOTT 

1. Isotropic eddy viscosity models 

2. Reynolds stress & scalar flux algebraic equation models 

3. Second moment transport equation models 

4. Multiple- scale models for compressible turbulent flows 

5. Bypass transition models 

6. PDF models for turbulent reacting flows 


PROGRAM SUMMARY 
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Isotropic eddy viscosity models 


• Objective 

0 To examine the deficiencies of existing models 
To develop better eddy viscosity models 

• Current status of existing k — e eddy viscosity models 

2 k 2 

-UiUj = u T {U itj + U jti ) - - k8ij , v T = Cpf^ — 

~ = T (fc) + - e + W.C., = t (c) + P< e > - D (e) + W.C. 

J-sL £/l 

They are not tensorially invariant due to f fl (y + ), W.C.(y + ) 

$ Model constants are not consistent for flows with and 
without wall 

Normal stresses may violate realizability 
^ Do not work very well for flows with pressure gradients 

• Development of a Galilean-, tensorially invariant, realizable, k — s model 

<0> New damping function f^k/Sv) is proposed to remove the depen- 
dence on y 

$ New dissipation £ equation is introduced to give better 
response to pressure gradients 

0 Consistent model coefficients for all flows 

Realizability of the normal stresses is guaranteed 

Modified wall function for cases with pressure gradients 


7 



• CMOTT k — e eddy viscosity model 


2 k 2 

-muj = vr{Ui,j - Uj,i) - - k6ij , v T = 


Dk 

Dt 

Ds 

Dt 


= Tk + Pk - s 


^ = T. + Ci fiS e-C 2 


k 4- 




<> are functions of R = k/Su, which is tensorially invariant 

<> = Aq+a 1 t/ . which ensures realizability for normal stresses 

0 $ represents the effect of inhomogeneity 

k 2 k 4 

$ = Vfc Vfc + 62— V 5 Vfc + 63 ^-V 5 VS 

£ £ z 


• Validation 
Flows: 

0 Channel flows 

<£> Boundary layer flows with and without pressure gradients 
<C> Planar jet, round jet and mixing layer 
<> Backward-facing step flows 

^ Complex flows related to industrial applications 
Models: 

^ Launder- Sharma, Lam-Bremhorst, Chien, 
Nagano-Hishida, . . . 


k — u model (Wilcox) 
0 CMOTT k - e model 
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Algebraic Reynolds stress models 


• Objective 

0 To examine the deficiencies of existing ARS models 
0 To develop better ARS models 

• Current status of ARS models 

^ Second-order closure based ARS models (Rodi, 1980) 


UiUj 

~k~ 


( P 



—UiUkUj'k ~ UjUkUi,k ~ p{p,i u j V,j u i) 2v u i,k u j,k 


Comments: 

* Assumption: upuj/k = Const., (uiUjUk),k = ( kui ),i = 0 

* Numerical difficulties 

0 Pope’s explicit ARS model (2-D flows), Taulbee’s ARS model (3- 
D), Gatski and Speziale’s ARS model 

^ Other methods: RNG, DIA and invariant theory 
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General constitutive relations from invariant theory 


miLj = + ^ a 2~~{Uij + Uj ti - -Ui ti 6ij) + 2 + Uj ti - -Di^y) 

+ 2a&^(Ui t kUj'k ~ -n 2 5,_ 7 ) + 2 a 7 — (Uk,iUk,j - -n 2 $y) 

S o S o 


+ 2a$-^-(Ui t kUj k + Ulp s * — — Ds^) + 2o-\Q—^-{U k ^U k j + U k ,jU ki 
+ 2ai2-^-{Ul k U? k - + 2aiz — (U k i U k j - -Ht£iy) 

+ 2au^-(U i , k Ui, k U? ij + U j%k U hh Ul - |n»^) 

T ( 6 2 

+ 2a 16 + UuJJlMi - 3 n ^«) 

+ 2a^(u itk Ui. k vl m ui„ + u^ihtVlJj?^ - |n 7 « 0 ) 


• EDT and realizability constraints (Reynolds, Lmnley) 

• CMOTT algebraic Reynolds stress model 

mvj = hs^ - +ic£(-s: h ni j + n-^sy 

o S c 


where 


k t t + Ujkj — [{v + U{iLj U{j e. 

<r k 

^ s sP 

£,t + Uj£ t j = [(v + — )£ t j]j — C c i—UiUjUij — C c 2 


c. 


1 

A 0 + A* 


u-k ’ 


C 2 


\/l-9C2(^) 2 

C 0 + 6*1*0^ 


k 2 

u t = Cu — , A 0 = 6.5, C 0 - 1.0 
£ 

Cel = 1-44, C c2 = 1.92, <T k = 1, a c = 1.3 
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• Validation 


0 Rotating homogeneous shear flows 
0 Backward-facing step flows 
0 Confined jets 

0 Complex flows related to industrial applications 



Configuration of rotating homogeneous shear flow 

10 
8 
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St 

Evolution of turbulent kinetic energy with time. 

: present model; — : SKE; • : LES 
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SWIRLING FLOW 


LARGE EDOY SHEAR 

WAKE REGION BETWEEN JETS 
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DEVELOPED 
DUCT FLOW 
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- SKE 

- RRSAE 
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Scalar turbulence model 

• Objective 

To improve the predictive capability of current scalar 
turbulence (0 2 - £$) models 

0 A new scalar flux constitutive relation 

<£ A new scalar dissipation rate model equation 


UiO = - Ca — (-) 1>/2 0,t +^(-) 1/2 (a 2 Ui, j + a 3 Uj y i)&,j 

€ T £ V 


Ui Wj =!' ^-‘0 J ).J + C ‘> “ s + - c « T 


_ (2 + 2r + 0.5r 2 ) 

A 26 + 3.2?7 2 + 2£ 2 

St = V®7Q, l S= y /2S^S~, r] = Sk/e, ^ = f(i) 1/2 Sr, r = 

C'ei = Ci — 0.13, C 02 = 0.63, C $ 3 = Cj — 1, <7* = 1.0, = 1.8 
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Tavoularis & Corrsin Expt (1981) 




Flat plate boundary layer with constant surface temperature 



y 
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• Validation 

Homogeneous turbulence subjected to dO/dy 

<0> Homogeneous turbulence subjected to dUfdy , dQ/dy 

0 Flat plate boundary layer with constant surface 
temperature 

• Work in progress 

^ Model assessment for different scalar boundary conditions 
0 Model extension for integration to the wall 


Sirivat & Warhaft Expt (1983) 
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Second Order Closure Models 


= Ta + Pi, + n?.* pid + ng«»" - |rf y 

• Objective 

0 To assess existing models 

0 To find the direction of improving closure models 

• Basic model forms 

n§ ** pid = i), 

ng**”" = F„(v-v-,v,k,e), 

Tij Fij (("Ujlij ) t A: 5 c) 

• General comments on second order closures: 

0 The model, I3^ apid , is relatively well developed compared with 
other terms 

0 The model, H^ eturn , is least developed 

0 A Galilean and tensorially invariant second order closure model 
has not been well developed yet 

§ All models have large errors near the wall, especially in the buffer 
layer; therefore, for engineering application, the wall function ap- 

proach is suggested at the present time 
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• Application of realizability to IP and LRR models 



m in 
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Multiple scale k-e model 


• Objective: 

^ To consider the effect of a non- equilibrium energy 

spectrum on eddy viscosity for compressible turbulence 

• Approach: 

0 Use multiple scale concept introduced by 


□ Large-Scale 


-D k p d PT,dk p . ,du, ? 

P Dt dy^ + a~ ) dy ] + ~ pCp + fCl 


De p _d. u T de p e p ,du. 2 _e „ 2 

p Dt - dy l(p+ Ct ’T r f ‘ T{ a^ ) ~ Cp2 % + fC2 

• fcj - exchanges between the turbulent kinetic energy and internal energy 

• fc 2 - increased spectral energy transfer due to compressibility effects 


□ Small Scale 


-D k t d Mr \ dk t , — 

”W~ a? K ' ,+ +Ph ~ **• 

K t 

_D e t d u T de t _e t e p _e t 

Pj= 7 = M + — )-x- + Ctxp-sf- - Ct 2 p T 

D t dy a~ dy k t k 


□ Eddy Viscosity 


y T m pul « p(k p + k t \ 2 


i (k p + k t ) 2 
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Model Evaluation 


• Turbulent Shear Flow 



TURBULENT FREE SHEAR LAYER TURBULENT BOUNOARY LAYER 


. Shock/Turbulent-Boundary- Layer Interactions 
0 transonic flow 



BUMP FLOW BUMP-IN-CHANNEL BUMP-IN-CHANNEL 


0 supersonic flow 



SHOCK REFLECTION RAMP FLOW 


Compressible Turbulent Shear Flow 



TURBULENT FREE SHEAR LAYER TURBULENT BOUNDARY LAYER 
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Bypass transition models 


• Objective: 

Develop transition models for flows with free stream 
turbulence 

• Approach: 

^ Using K-e model as the base model 

0 Introduce effective intermittency to either the eddy 
viscosity or the k-£ model equations 

Transitional Boundary --yer 



Transilionol Boundary Layer 


Tu-3% 
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PDF modeling of turbulent reacting flows 


• Objective: 

0 Develop models that can accurately simulate finite chemical reac- 
tions in turbulent flows. 

^ Develop and validate independent PDF models. 

0 Technology transfer. 

• Approach: 

^ Joint pdf for scalar compositions. 

^ Moment closure schemes for velocity field. 

0 Develop hybrid solver consisting of Monte Carlo method and 
finite-difference/ finite- volume method. 


BURNER TET 




Without PDF Method 



With PDF Method. (CfflOTT) 


^ _ 



* 
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